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Abstract 

We find that fermion determinants are not generally positive in a recent class 
of constructions with explicit lattice supersymmetry. These involve an orbifold 
of super symmetric matrix models, and have as their target (continuum) theory 
(2,2) 2-dimensional super- Yang-Mills. The fermion determinant is shown to be 
identically zero for all boson configurations due to the existence of a zeromode 
fermion inherited from the "mother theory." Once this eigenvalue is factored 
out, the fermion determinant generically has arbitrary complex phase. We 
discuss the implications of this result for simulation of the models. 
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Introductory remarks. Models with exact lattice supersymmetry have been dis- 
cussed in the literature by a few groups. For example, latticizations of super- Yang- 
Mills mill, supersymmetric quantum mechanics 011] j the 2d Wess-Zumino model 
E] , and direct constructions in the spirit of the Ginsparg- Wilson relation — as sug- 
gested by Liischer jU] and worked out in noninteracting examples [HIH] — have all been 
considered.^ In this letter we will be interested in the super- Yang-Mills constructions 
that lead to a Euclidean lattice theory PP. The method of building such models is 
based on deconstruction of extra dimensions jTUl The corresponding interpre- 
tation in terms of the world-volume theory of D-branes has led to the latticizations 
of 2d, 3d and 4d supersymmetric gauge theories. These lattice constructions are all 
arrived at by orbifold projections of supersymmetric matrix models; i.e., in each case 
we quotient a matrix model by some discrete symmetry group of the theory. Degrees 
of freedom that are not invariant with respect to the combined action of the orbifold 
generators are projected out.^ Thus, throughout this letter we will have occasion to 
speak of "orbifolded" matrix models and "nonorbifolded" matrix models. 

A major motivation for efforts to latticize supersymmetric models is that some 
nonperturbative aspects of supersymmetric field theories are not accessible by the 
usual techniques, such as holomorphy. One hope of a lattice supersymmetry program 
of research is that it would lead to, e.g., simulations that would provide further data 
on supersymmetric field theories, especially those that include super- Yang-Mills.^ 
With exact lattice supersymmetry the target (continuum) theory may be obtained 
in a more controlled fashion. Indeed, in some cases it may be obtained without the 
need for fine-tuning [U |2j . 

In this enterprise, it is of great practical importance that the fermion determinant, 
obtained integrating over the fermion degrees of freedom in the partition function, be 
positive. For let be the lattice bosons in the theory, and the lattice fermions. 

^An approach very similar to [H| has been apphed in [H], yielding slightly different expressions. 
^For a detailed discussion, we refer the reader to p. 

■^For a recent review of existing work on this broad topic, and a complete list of relevant references, 
see [E]. 
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A positive^ det M(0) for all allows us to unambiguously calculate in an equivalent 
bosonic theory with action Ses{4>)- Expectation values of operators 0{(j)) may be 
obtained using this action: 



Fermionic correlators will simply involve M~^(0) in the operator of interest. In the 
case of positive detM(0), the techniques of estimation by Monte Carlo simulation 
are robust; otherwise the situation is murkier. 

We now summarize the content of our work: 

• In this letter we show that the lattice theory with (2,2) 2d super- Yang-Mills as 
its target P, obtained from orbifolded supersymmetric matrix models, possesses 
a problematic fermion determinant. Due to a zeromode fermion, det = 0, 
i.e., for all boson configurations. 

• We then suggest how the zero eigenvalue can be factored out in a controlled way 
in order to exhibit the determinant for the other fermions. For the orbifolded 
matrix models, we carry out this factorization (numerically) in the special case 
of a U{2) gauge theory, and show that it is robust. As an example, we present 
our results for a 2 x 2 lattice. 

^Here we assume ordinary finite dimensional Berczin integration for the fermion coordinates; thus 
we avoid subtleties of the Weyl determinant present in continuum theories . The class of 

theories studed here result from a 4d ^ Od dimensional reduction, where the Pfaffian one obtains 
has an equivalent determinant form I15j . 

^Strictly speaking, we require only positive semi-dcfinitc det A/(0), provided detA/((/)) — oc- 
curs only for some subset of the boson configurations; this subset has measure zero. However, if 
detM((/)) = 0, i.e., for all boson configurations, then Z is not well-defined. This latter situation is 
what occurs in the theory that we study here. 





(3) 
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• We further validate our method by applying it to nonorbifolded U{k) supersym- 
metric matrix models, which also contain ever-present zeromode fermions, cor- 
responding to the U{1) "gaugino" in the decomposition U{k) D f/(l) x SU{k). 
We correctly reproduce the fermion determinant for the nonorbifolded SU{k) 
matrix models by our factorization method. As an example, we present our 
results for f/(5) D f/(l) x SU{5). 

• Once the zeromode fermion has been factored out of the orbifolded matrix mod- 
els studied here, we find that the remaining product of eigenvalues is generically 
nonzero with arbitrary complex phase. 

• We explain how this is not in conflict with results in nonorbifolded SU{k) 
supersymmetric matrix models. (For example, in an appendix of |TH] it was 
shown that the product of eigenvalues for the fermion matrix is positive semi- 
definite.) It is found that the orbifold procedure lies at the heart of this matter. 

• We conclude with a discussion of the implications of out results for lattice 
simulations of the latticized (2,2) 2d super- Yang-Mills theories. 

Zeromode fermion. Here we will focus on the orbifolded supersymmetric ma- 
trix models that have as their target theory (2,2) 2d super- Yang-Mills with U{k) 
gauge group. The "mother theory" is a nonorbifolded U{kN'^) supersymmetric ma- 
trix model. The "daughter theory" is obtained by orbifolding the "mother theory" by 
a Zjv X symmetry group, leaving intact (among other things) a U{k)^^ symmetry 
group that will become the gauge symmetry of the N x N lattice theory. The lat- 
tice theory is obtained by studying the "daughter theory" about a particular boson 
configuration that is a stationary point of the action; i.e., it is a point in the moduli 
space of the "daughter theory." 

The orbifolded matrix model discussed here has been described in detail in [1 ; we 
refer the reader there for further details. For our purpose it suffices to note that in 
the U{k) case the theory contains^ bosons Xm, Vm that are k x k complex matrices. 
The bosons x^, Vm may written in terms of a Hermitian basis 

T^G |y|lfc,T'^|, (T'^)t = T^ x^ = xii,T^ = (4) 

^Our notation is standard: m labels sites on a 2d square lattice, so that m G Z x Z, with Z the 
set of integers. 
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It is always possible to choose the such that 

Tr {Tf'T'') = 26^". (5) 

Furthermore we define 

Tr {T^'T'^T") = 2^1 1^""" (6) 
and note that (underlining implies all permutations are to be taken): 

One finds that the fermionic part of action is given by 

Sf = P^) ■ ■ (^^] . (8) 



) /7~V m'r^m/ m,n 1 /-p 

g^vk V^n 

Here a^, A^, are the lattice fermions, with upper index corresponding to the 
basis T'* introduced above. The fermion matrix M^^^ is given by (sum over u implied 
in the entries, i,j unit vectors): 



m,n 



m,n m 


m,n— ? m 


^ff^l^P 7,1' J 
m,nym+i 




j-Hup jjU _ 
''m,ny m 


_ fi^P^ TrjV 
m,n— jyin 


fpi^p _ 





(9) 



Here we have introduced the compact notation 

C"n = 5m,„i'^'^''. (10) 

The fermion zero mode is easily established. We consider fermions of the form 

Vn. (11) 



r 



Then in this case 

E—iy (^tXVp _ .p.pl' \ \p _ \ ^ —ly /y./ivO _ ^P.0U N ^ 

n n 

= <a5^(5„,„-5„,„_0 = 0. (12) 

n 

By similar arguments 

E ^m(Cfn - Cn-j) = 0. (13) 
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Thus (lllj) is an eigenvector of M with eigenvalue zero. It follows that det M (x, y) = 0, 
for all configurations of bosons. This clearly poses a difficulty for defining Scs{x,y)- 
We will shortly address a method to factor out this zero eigenvalue, and study it in 
some detail for the U{2) case. However, we first make a few remarks on the existence 
of the indentically zero eigenvalue. 

The zeromode fermion is nothing but the "zero-momentum" mode of the Fourier 
transform of a given A^: 

^-^o = ^E^- (14) 

m 

In the formalism introduced in [Ij, appears in the superfield 

Here d!^ is an auxiliary boson and 9 is an odd (Grassman) superspace coordinate. 
The zero-momentum part of this supermultiplet is just 

A = 5^A^ = A + zrf^, d^j^J2< (16) 

m m 

That is, the zeromode fermion is in a multiplet that contains just itself and an auxil- 
iary boson. It is a lattice version of a Fermi multiplet jT7|. Thus there is no physical 
zeromode boson that corresponds to the zeromode fermion. 

The existence of the zeromode fermion can be understood in terms of the "mother 
theory;" i.e., the nonorbifolded matrix model. There, the bosons x,y and their con- 
jugates are understood in terms of a "vector boson" v. 



V = VmCTm = Vq + iv CT = ^2 [ _ (17) 




where Vm are Hermitian matrices that are Lie algebra valued in U{kN'^): 



\^m)iiii/,jfM'u' 



v^{T%,,,^,,, t,j = l,...,k, fx,u,fx',u' = 1,...N. (18) 



Note that the indices of generators have been written in such a way that the U{k)^^ 
subgroup has been manifestly factored out. One should think of one U{k) factor 
"living" at each site. Similarly, the fermions of the mother theory are given by 

i^i.^.^'u' = ^"(T")w,.^v', r = QI) (19) 
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with a corresponding expression for ijj. The fermion action in the mother theory takes 
the form 

Sf = -^Tt [ipa^[vm,i']) ■ (20) 
Now note that the diagonal f/(l)diag C U{kN'^) fermions do not appear Sp- That 

is, 

9 (21) 

and ip^ disappears because of the commutator. Since ip^ is a two component fermion, 
it gives two zeromode fermions of the mother theory, independent of the boson con- 
figuration Vm- The boson action takes the form 

Sb = {[Vm,Vn][Vm,Vn]) ■ (22) 

Here again, the U (l)diag boson disappears from the action because of the commutators. 
It is a zeromode for all boson configurations. 

Most of the fields of the mother theory are projected out in the orbifold construc- 
tion. The projections depend on charges with respect to a f/(l)ri x f/(l)r2 global 
symmetry group. As it turns out, the A" fermions are f/(l)ri x U{l)r2 neutral. It fol- 
lows that only the diagonal parts with respect to the "site" indices (/iz/, /xV) survive 
in the orbifolded matrix model ("daughter theory"): 

(T )i^u,j^j,'u' A (T )i^u,j^j,u ■ (23) 

But this includes the zeromode fermion of the mother theory. On the other hand, 
none of the components of v are f/(l)ri x U{l)r2 neutral. It follows that, in the orb- 
ifolded matrix model, they are all off-diagonaF with respect to the indices (/xi/, /iV). 
Hence the zeromode bosons of the mother theory are projected out in the orbifolded 
theory. Similarly, none of the ^" fermions of the mother theory are ?7(l)ri x U{l)r2 
neutral, so that the second zeromode fermion C,^ of the mother theory disappears in 
the projection. 

Deformed U{k). The zeromode eigenvalue of the orbifolded matrix model can 
be factored out as follows. We deform the fermion matrix according to 

M ^ M + elN^ (24) 

^In the lattice theory, these are interpreted as hnk fields beginning at the site labeled by {ii, v) 
and ending at the site labeled by (/^', v'). 
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where Nf is the dimensionahty of the fermion matrix and e <^ 1 is a deformation 
parameter that we will eventually take to zero. We factor out the zero mode through 
the definition 



If this deformation is added to the action, it explicitly breaks the exact lattice 
supersymmetry and gauge invariance. This infrared regulator could be removed in 
the continuum limit, say, by taking ea <^ A^~^. Noting that L = Na is the physical 
size of the lattice, the equivalent requirement is that e ^ L^^ be maintained as a — > 0, 
for fixed L. Thus in the thermodynamic limit {L oo), the deformation is removed. 
The parameter e is a soft infrared regulating mass, and is quite analogous to the soft 
mass n introduced by Cohen et al. in their Eq. (1.2) to control the bosonic zeromode 
of the theory. In the same way that a finite n does not modify the final results of the 
renormalization analysis of Section 5 of our e does not modify the result of the 
quantum continuum limit. The essence of the argument is that we have introduced a 
vertex that will be proportional to the dimensionless quantity ^ffea^ ^ gla?' / where 
g2 is the 2d coupling constant. Such contributions to the operator coefficients Co in 
Eq. (5.2) of Ij vanish in the thermodynamic limit and because the target theory is 
super-renormalizable, we are assured that the perturbative power counting arguments 
are reliable and the correct continuum limit is obtained. 

Of course it remains to study the convergence of detM(e) — > detM(O). Indeed, 
for U{2) lattice theory (i.e., the orbifolded matrix model with f/(2)^^ symmetry) we 
have performed this analysis numerically for a large number of boson configurations 
drawn randomly from a Gaussian distribution, as will be detailed below. We find that 
the convergence is rapid and that a reliable estimate for det M(0) can be obtained 
in this way. Furthermore, we find that for e <^ 1, the phase of det M(e) quickly 
converges to a constant value, and that it is uniformly distributed throughout the 
interval (— tt, tt] for the random Gaussian boson variables. 

As a check of our method, we have studied also the analogous deformation of the 
nonorbifolded U{k) supersymmetric matrix models, where the two zeromode fermions 
in tp^ are present. Indeed, as will be discussed below, we find that 



M(0) = limM(e), M(e) = e-^'^f{M + eljy,), 
^ detM(O) = lime"Met(M + eljvJ. 



(25) 



lime ^ det(M[/(fc) + elNj) = det Msu{k)- 



(26) 
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That is, we obtain the determinant of the nonorbifolded SU{k) supersymmetric ma- 
trix model, which has the zeromode fermions factored out. With appropriate conven- 
tions, det Msu{k) is positive semi-definite.^ 

Deformed U{2). Here we speciahze to the orbifolded matrix model with U{2)^^ 
symmetry, which becomes the U (2) gauge invariance of an iV x iV lattice theory. In 
this case we take 

T^e{l2,an, x^ = xl, + xl,a\ |/m = Z/^ + (27) 
Then the fermion matrix is given in Q, where in the present case 

^000 ^ ^oOO ^ ^gm ^ ^ab ^ ^abc ^ -^abc^ ^28) 

The lattice theory is obtained by expansion about a point in moduli space: 

•^m ~ 7^ + ■ ■ ■ ) Vm ~ 7^ + ■ ■ ■ 5 (29) 
av2 ay/ 2 

where ■ ■ ■ represent the quantum fluctuations. For this reason, in our study of 
det M(e) we scan over a Gaussian distribution where x'^,y'^ have a a nonzero mean 
l/ay/2 = 1. The remainder of the bosons are drawn with mean zero. All bosons are 
taken from distributions with unit variance. 

In Figs. HI and El we display ln|detM(e)| and argdetM(e) versus e for the case 
of = 2, the smallest lattice possible, of size 2x2. Each line corresponds to a 
different random draw. It can be seen that the convergence is rapid and reliable. As 
a check, we have computed the eigenvalues of the undeformed matrix M, using the 
math package Maple, for the same set of random boson configurations. We find that 
in each case the product of nonzero eigenvalues agrees with det M(0) in magnitude 
and phase to within at least 5 significant digits. This conclusively demonstrates that 
the complex determinant is not an artifact of the deformation, but is a property of 
the nonzero eigenvalues of the undeformed matrix M. 

For a set of 10^ draws on the bosons of this 2x2 lattice, we have extrapolated to 
e ^ and binned argdetM(e) over its range, with bins of size vr/lOO. In Fig. 01 we 
show the frequency for each bin, as a fraction of the total number of draws. In the 

^For some conventions on the SU (fc) generators and the overall coefficient in Sp of Eq. (|20|l it is 
possible that a constant phase (i.e., independent of boson configurations) may be present. However, 
this factors out of Q and is of no concern. 
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Figure 1: ln|detM(e)| versus e for a sequence of random draws. These results are 
for the U{2) lattice theory, with 2x2 lattice. 



extrapolation, we decreased e by powers of 10 until the change in both In | det M(e)| 
and argdet M(e) were both less than 10~^ per decade. This was done for each draw 
to get a reliable estimate for det M(0). In each of the 10^ draws the 10"'^ per decade 
criterion was reached well before e = 10^^°. 

To summarize, once the zeromode eigenvalue is factored out, the product of the 
nonzero eigenvalues has arbitrary phase. Consequently an ambiguity exists in defining 
Scs{x, y) for the orbifolded matrix model theory. We are presently exploring whether 
or not this can be overcome for the purposes of simulation. At present, however, all 
we can say is that this difficulty is rather troubling. 

Comparison to nonorbifolded U{k) and SU{k) matrix models. As mentioned 
above, the nonorbifolded supersymmetric U{k) matrix models contain two ever- 
present zeromode fermions in the ip^ that appears in ()21|) . Then detM(7(fc) = 0. 
However, it is easy enough to just work with the nonorbifolded SU {k) matrix model, 
so that ip^^ip^ are never in the theory to begin with. Then with appropriate con- 
ventions det Msu{k) > 0. A proof of this result has been given in an appendix of 

As a test of our method, we have verified numerically, for a sequence of 
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IQ-iO lO^S 10^6 10^^ 10^2 iqO 

e 



Figure 2: (j) = argdctM(e) versus e for a sequence of random draws. Note that 
crossing over the boundary of the domain (— vr, tt] to an equivalent point within that 
domain is indicated by the nearly vertical lines. These results are for the U (2) lattice 
theory, with 2x2 lattice. 
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Figure 3: Frequency distribution F{(j)) for = argdetM(O), for 10^ random (Gaus- 
sian) draws, binned into intervals of vr/lOO. The distribution of is seen to be nearly 
uniform. These results are for the U{2) lattice theory, with 2x2 lattice. 

random Gaussian draws on the bosons v"^ that appear in ((20)). In Fig. HI we show the 
quantity 



A(e) = In [e"2| dei{Mu{^k) + eljV;)|] - In | det M, 



'-SU{k)\ 



(30) 



as a function of e for the case of = 5. Indeed it can be seen that the convergence as 
e ^ is quite rapid. We find that |A(10-3)| < 10"^ and |A(10-^)| < lO'^ as a rule. 
In Fig. El we show the quantity 

0(e) = arg [e'^ det{Muik) + eljv,)] (31) 

as a function of e, again for k = 5. The value for SU (5) in our conventions is = 0. 
It can be seen that a rapid convergence to this value is obtained. 

One might wonder how the generic phases of Fig. El are obtained in the orbifolded 
matrix models, given the positivity of the fermion determinant in the nonorbifolded 
SU{k) matrix models. Firstly, one should note that the proof given in ^B] shows only 
that for each eigenvalue e of Msu(k) there exists also an eigenvalue e*, and that these 
always come in pairs. The proof relies essentially on the relation 



cT"(w^ + iv"-<T)a^ 



(t;°)*. 



(32) 
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Figure 4: A(e) versus e for a sequence of random draws in the nonorbifolded U{5) 
matrix model. 
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Figure 5: (p versus e for sequence of random draws in the nonorbifolded U{5) matrix 
model. 
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On the other hand, the x orbifold action on the mother theory bosons is 
generated by 

)ipfi,jp'p,' > e ^ V e ^ VtpyiT )ip,ujfi'v'^y' pi (33) 

where VL = diag {lj, u^, . . . , u^), u = exp{2ni/N). Thus because of the exp(±i7rcr^/A^) 
factors, the orbifold projection does not commute with the operations of the proof 
given in JE]; i.e., Eq. ((221) • The fermion matrix of the projected theory — i.e., the orb- 
ifolded matrix model — lacks many of the eigenvalues of the mother theory; it should 
come as no surprise that not all eigenvalues are removed in pairs (e, e*). After all, 
we already know that only one of the zero eigenvalues is removed from the U{kN'^) 
mother theory. For this reason it is not at all contradictory that the projected theory 
has a product of nonzero eigenvalues that is not positive, nor real. 

Discussion. The existence of a zeromode fermion in the orbifolded matrix mod- 
els has been reliably handled by our deformation method, and the approach is easily 
implemented numerically. Analytic methods in nonorbifolded supersymmetric matrix 
models have also involved deformations of the theory to render quantities of interest 
well-defined ^S|; however, these involve the introduction of auxiliary fields, which 
are expensive to implement in a simulation. For this reason, we prefer the method 
described here. In spite of the complex fermion determinant in the orbifolded matrix 
models, our approach allows for a detailed study of the partition function by Monte 
Carlo methods, analogous to what has been performed in J3] for nonorbifolded su- 
persymmetric matrix models. It would be interesting to see how expectations based 
on continuum results might be realized in the present context. For example, in ^3] it 
was shown that 4d A/" = 1 pure super- Yang-Mills has a positive fermion determinant. 
Since the (2,2) 2d target theory studied here is a dimensional reduction of this model, 
it is rather surprising that we find arbitrary complex phase. However, it remains to 
be seen whether or not the lattice action described here possesses reflection positiv- 
ity, or whether this is only obtained in the continuum limit. Study of this issue is in 
progress. 

We suspect that the difficulties faced here may be, broadly speaking, related to 
those faced in defining the phase of the fermion measure in the attempts to realize 
chiral gauge theories on the lattice with an exact chiral gauge symmetry ^01 IIHI 
1^ . In the present context, we have an exact chiral fermionic symmetry: lattice 
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supersymmetry acts on the Ain,^ni but not on the a^, Pm that appear in (jH)). A 
resolution of one problem may lead to answers for the other. 

In our opinion, the deconstruction approach to lattice supersymmetry remains 
an exciting topic, whatever difficulties may face attempts to simulate the theory. 
For example, the relatively simple systems described here provide another context in 
which the difficulties^ that plague the complex Langevin approach fI^ 1^ and other 
complex action techniques might be studied. 

We are presently exploring whether or not the complex phase can be overcome for 
the purposes of simulation. A typical approach would be to compute averages of an 
operator O from the re-weighting identity: 

iO) = (34) 

Here, = argdetM(O), as above, and "p.q." indicates phase- quenching: expectation 
values are computed with the replacement det M(0) — > | det M(0)|. Thus the effective 
bosonic action 

5'^g = 55-ln|detM(0)| (35) 

is used to generate the phase-quenched ensemble by standard Monte Carlo tech- 
niques. However, it is well-known that this suffers efficiency problems: the number of 
configurations required to get an accurate estimate for, say, (exp(20))p^ , grows like 
exp(AF ■ N'j). Here, AF is the difference in free energy densities between the full 
ensemble and the phase-quenched ensemble. Recall that Nj is the dimensionality of 
the fermion matrix. On the other hand, if the phase-quenched distribution is highly 
concentrated near one value of 0, the phase essentially factors out of the partition 
function and efficient simulations can be done using the phase-quenched ensemble. 
Thus it is of interest to study the distribution of in the phase-quenched ensem- 
ble rather than the Gaussian distribution used here. Research in this direction is 
currently in progress and we hope to report on it soon. 
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